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SUPERSONIC NATE BRAG OF NONLIFTING SWEPTBACK 
TAPERED WINGS WITS MACH LINES BEHIND 
THE LINE OF MAXIMUM THICKNESS 
By Kenneth Margolis 


SUMMARY 

A theoretical investigation of the supersonic wave drag of nonlifting 
sweptback tapered, wings having thin symmetrical double-wedge airfoil 
sections with maximum thickness at 50 percent chord has teen presented in 
NACA TN No. l4H8. The present paper extends the investigation to include 
"supersonic" maximum-thickness edges ; that is, the flight velocity component 
normal to the line of max im um, thickness is supersonic. This condition 
exists at Mach numbers for which the' Mach lines have angles of sweep 
greater than that of the line of maximum thickness. 

For wings of equal root tending stress (and hence of d i fferent aspect 
ratio) and given sweepback, taper increases the wing wave-drag coefficient 
at Mach numbers for which the maximum-thickness line is moderately supersonic 
and has negligible effect at higher Mach numbers. This trend is similar to 
that evidenced by the effect of high aspect ratio for given sweepback and 
taper ratio. Comparisons on the basis of constant aspect ratio for given 
sweepback, however, indicate a decrease of the drag coefficient with taper 
at Mach numbers corresponding to moderately supersonic maxi mum— thickness 
lines and a negligible effect due to taper at the higher Mach numbers. 

For given taper ratio and aspect ratio, increased sveepbaCk increases the 
drag coefficient at Mach numbers for which the maximum-thickness line iB 
supersonic. 


INTRODUCTION 

In reference 1, equations are derived and calculations are presented 
for the supersonic wave drag of sweptback tapered wings with thin symmet- 
rical double-wedge sections with maximum thickness, at 5 0 percent chord 
(that is, rhombic profiles) at zero lift. The range of Mach number con- 
sidered in reference 1 is between 1.0 and the value corresponding to the 
condition where the Mach lines are parallel to the maximum-thickness line, 
that is, the Mach number range corresponding to "subsonic" maximum- 
thickness lines. (When the flight velocity component normal to an edge is 
subsonic, the edge is termed a "subsonic" edge. J This condition exists 
at Mach numbers for which the Mach lines have angles of sweep leBS than 
that of the line of maximum thickness. 

The present paper completes the wave-drag investigation of nonlifting 
tapered wings by extend i ng the Mach number range to include "supersonic" 
maximumr-thickness lines; that is, the region where the Mach lines have 
angles of sweep greater than that of the line of maximum thickness. 
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Accordingly, equations are derived for the supersonic wave drag of 
sweptback tapered wings with thin rhombic sections at zero lift. Calcu- 
lations are presented for some representative plan forms. As in 
reference 1, the airfoil sections and wing tips are chosen parallel to 
the direction of flight, and the angle of sweepback is referred to that 
of the line of maximum thickness. For purposes of— completeness, the 
results obtained in reference 1 for the lower supersonic Mach numbers are 
included in the discussion and conclusions of the present investigation. 


SYMBOLS 


x, y, z 

V 

P 

Ap 

4 

M 

p =\/m® 

dz/dic ' 

a 

c 

t 

A 

“0 

m l 

Dig 

b 


Cartesian coordinates 
velooity in flight direction 
density of air 
pressure increment 
dynamic pressure 
disturbance-velocity potential 
Mach number 



slope of airfoil surface, measured in flight direction 

root semichord, measured in flight direction 

chord length at spanwise station y, measured in flight 
direction 

maximum thickness of section at spanwise station y 
angle of sweep of maximum- thickness line, degrees 
slope of line of maximum thickness (cot A) 
slope of wing leading edge 

/ mj_m 

slope of wing trailing edge f 

\2mi - 

span of wing 
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mo 

S 

A 

X 

°<3oo 

C *tip 


c d 

% 

C DtIp 

°D 


wing area 

aspect ratio <?» 

taper ratio (ratio of tip chord to root chord) 

section ware-drag coefficient at spanvise station y 
exclusive of tip effect 

increment in section wave-drag coefficient at spanwise 
station y due to tip 

section wave-drag coefficient at spanwise station y 

(°a» + 

wing wave-drag coefficient exclusive of tip effect 
increment in wing wave -drag coefficient due to tip 
wing wave-drag coefficient + Op^N 


Subscript s refers to conditions at root 


ANALYSIS 


The analysis is essentially that used in references 1 and 2. A 
brief outline of the basic equations is included for convenience. 

The assumptions of small disturbances and a constancy of sonic 
velocity throughout the fluid lead to the linearized equation for the 
disturbance-velocity potential cp 



^xx + ‘•Pyy + ?zz 


0 


(1) 


where M is the Mach number of the flow and the derivatives are taken 
with x-espect to the variables x, j , and z of the Cartesian coordinate 
system. On the basis of this linear theory a solution for a uniform 
semi -inf ini te sweptback line of sources is derived in reference 2. The 
pressure field assoc'<?.ted with this solution corresponds to that over an 
airfoil of wedge section. The pressure coefficient Ap/q. at a spanwise 

station y and point x along the wedge is. for P , 

m l 



k 


' HACA TN Ho 


1672 


Ap 

q. 


2 dz 
it dx 


m l 



•, x - mqP^y 

cosh. — 

'■"? 7 ■ rn-jx 


and, for 



Ap _ 2 dz mi _i x ~ y 

1 ^ 003 


(2a) 


(2b) 


where is the slope of -the leading edge of the wing, dz/dx is the 

tangent of the half^wedge angle (approx, equal to half,- wedge a n gl e since 


( 


the angle is small), p = \Jb £ - 1, and the origin of-Jdi© line source is 
taken at (0,0). In the region between the leading edge and the Mach cone 

1 ^ ^ \ i x - m-,P 2 y 

: is, •s-x h- y = m-ix) , the real part of cos" — , — -2: 

P ' p I y - m-jx I 

and equal to it. Equation (2b) then reduces to 


is constant 


Ap _ 3 dz m I 
q dx ; 2~o 7 

/ m l P ' 1 


(2c) 


The distribution of pressure over sweptback wings of desired plan 
form. and profile is obtained by superpositioh of solutions for wedge-type 
airfoils. Reference 1 adequately describes the superpositions necessary 
to satisfy the boundary conditions over the surface of a tapered wing of 
rhombic section, and therefore the procedure will not be restated herein. 
Figure 1 shows the distributions of sinks and sources for a tapered wing 
and identifies the system of axes and the symbols associated with the 
derivation of the drag equations. 

The disturbances caused by the elementary line sources and sinks 
are limited to. -the regions enclosed by their Mach cones. Figure 2 shows 
the Mach line configuration for the tapered-wing plan forms and indicates 
the regions of the wing affected by each line source, and sink. For 
purposes of simplif ication the tapered wings considered were restricted 
to those with no tip effects other than the effects each tip exerts on its 
own half of the wing. For a wing of taper ratio 0, ho tip effects need 
be considered since, the. Mach lines originating at the tip do not enclose 
any part of. the wing. - - - 
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The pressure coefficients obtained from superposition of solutions 
given in equations (2) are converted into drag coefficients by the 
f ollowing relations : 


C D 





1T.E. 


Ap 

q. 




L.E. 


(3) 


where b is the wing span, S is the wing area, dz/dx is the slope of 
the airfoil surface, and L.E. and T.E. denote leading edge and trailing 
edge, respectively. 


DERIVATION OF GENERALIZED EQUATIONS 


By appropriate superposition of solutions for wedge-type airfoils, 
the pressure field is obtained for a tapered wing with leading edge, 
trailing edge, and line of maximum thickness sweptback. The drag equations 
are derived for half of the wing since the drag is distributed symmet- 
rically over both halves. The induced effects of the opposite half -wing 
are represented by the conjugate terms in the integrands of the drag 
integrals . 

For a rhombic profile. 


dz t 
dx “ c 

where t/c is the section thickness ratio. The generalized equations for 
supersonic wave drag exclusive of tip effects are obtained as follows 
(see fig. 3 for information pertinent to integration limits) : 
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For p > — , 
m l 

ltSC Ib 3 = 
8(t/c) 2 
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where 


1 x.+ a + miJ3 2 y x +. a - miP 2 y 

A=0 ° ah ' ?i7T^TTTi + °° 8h ‘ J\T- m l(l+ a) 


1 x + atf 2 ? 1 X - HqP 2 ? 

B = cos r + cos — * 

P y + fflQX P | y: - e^qX j 


and 


Q 2 p 

x + a + m x p y -1 x + a - m^y 

C = 003 P |y + mi(x + a) | + 003 p |y - mi (x + a) | 


It should he noted that equations (4) give the drag for plan-form 
configurations where the tip is placed farther spanwise than the points 
of Intersection between the Mach lines and the trailing edge . (See fig. 3 • ) 
Deletion of certain integrals and appropriate changes in the y limits 
of other integrals may he made for configurations where the tip is placed 
nearer the root^ chord. These equations are evaluated and the resulting 
section wave-drag and wing wave-drag formulas for all tapered plan forms 
are presented in appendix A. 

As stated previously, the tapered wings considered have no tip effects 
other than those that each tip exerts on its own half of the wingj that is , 
the Mach lines from one tip do not enclose any part of the opposite half- 
wing. This condition is expressed mathematically as follows: 


Aspect ratio = 


2dmg 


4mp 


a(l + 1) (1 + X) (1 + m.±$) 


(5) 


/Tip chord \ 

where X is the taper ratio' ). It can he seen from 

\Eoot chord/ 

equation ( 5 ) that this simplification does not materially limit the ranfp 
of Mach number that may he considered. In fact, since these equations 
involve only Mach numbers corresponding to supersonic maximum— thickness 
lines, the limiting effect imposed in equation ( 5 ) is found to he negligible. 
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The wave-drag contribution of the tip is (fig. 3 ): 
For — < p < -L, 

TO/-. — — -r* ' 
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y+nigS' 


mg 


D dx dy 


l IL. 


m 1 (l+m 0 P ) mo 


(l+m^p ) -m^a y+m^a 


m i(i +E1 oP ) 


/]' n, 
fflg 


m^n^dCl+ii^p ) -201^^ Jm^Cl+n^p ) -m^a+py) 
nrpCl+mgP ) m.j_ 


D dx dy 


2 m„ 


) m o d - 


\/™o 2 P 2 - 1 


l+mgP 
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( 6 ) 
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where 


P P 

m 1 (x + a) - hIqI - p (y 


D = cosh 


-1 


m o d) 


m lP| y - m i(x + a)| 


and 


E = cos' 


x - d - m 0 P 2 (y - n^d) 

3jy - mox| 


For p > — , equation (6) is still- valid if /l - m-j_ 2 p 2 is changed to 

i/m^p 2 ~ and the inverse hyperbolic function cosh - ' 1 is changed to 

the inverse cosine function cos“ d . - 

Equation (6) was solved for section wave drag and wing wave drag, 
and the results are presented -in appendix B. The total wave-drag coef- 
ficients are then obtained by the following relations: 


G d = c d M + c d tip 

C D = C Doo + C Dtip 


(7) 


The value of 


Cjj^_ is found to be identically equal to zero for all 


cases satisfying the limitations imposed in equation ( 5 ) and, hence, 
Cp = Cp, for the tapered wings, considered. 


RESULTS AND DISCUSSION 


The formulas presented in the appendixes apply to all conventional 
taper ratios (0 ^ 1) . For the particular case of taper ratio 1 

(untapered- plan form) the equations presented in appendix A of reference 3 
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are in a more convenient form for calculation purposes. Calculations 
presented for wrings with, subsonic maximum-thickness lines are based on 
the formulae of reference 1. 

Calculations are made for some typical tapered plan forms and also 
for a family of tapered wings considered in reference 1; members of a 
family are characterized by constant sweepback of the line of maximum 
thickness and a constant value of the parameter "moment of wing area about 
the root chord divided by tho product of the root chord and the square of 
the root thickness." This area— moment condition is intended to imply that, 
to a first approximation, the root bending stress is the same for all 
members of any family having the same thickness ratio. (See fig. 5 of 
reference 1 for further details . ) 

Section wave drag.— Spanwise distributions of section wave drag for 
wings of taper ratio 0, 0.5, and 1.0 are presented in figure k for a Madi 
number of 3 , aspect ratio 2, and sweepback of 60°. These representative 
spanwise drag distributions for wings with supersonic m,?xinn3,m— thickness 
lines differ markedly from those obtained at lower speeds where the wing 
is swept well behind the Mach lines. All sections have positive drag at 
speeds corresponding to supersonic max i m um-thickness lines, whereas, at 
the lower Mach numbers, outboard sections experience negative drag. 

(See figs. 6 to 10 of reference 1.) 

Wing wave drag .— Variations of wing wave— drag coefficient with Mach 
number for constant-aspect— ratio wings of taper ratios 0 and 1.0 and 6o° 
sweepback are presented in figures 5(a) and 5(b) for aspect ratios 2 and 3* 
respectively . For wings of constant aspect ratio, taper increases the 
drag coefficient at Mach numbers for which the max imum-thi ckne s s line is 
substantially subsonic, decreases the drag coefficient in the intermediate 
range, and has negligible effect when the maximum— thickness line is highly 
supersonic. For comparison purposes, the result for the straight wing of 
infinite aspect ratio (two-dimensional case) is included in figures 5 
to 7- A typical variation in wing wave-drag coefficient with Mach number 
is shown in figure 6 for wings of equal root bend in g stress. (The data for 
the curves between Mach numbers 1 and 2 are taken from fig. 12 of 
reference 1. ) Taper is shown to reduce the wing wave-drag coefficient at 
Mach numbers for which the maximum-thickness line is substantially subsonic, 
increase the drag coefficient in the intermediate range, and have little 
effect when' the maximunt-thickness line is highly supersonic. These trends 
are similar to the ones shown by the effect of high aspect ratio on the 
wave— drag coefficient of wings for a given taper ratio. (See fig. J.) 

It must he remembered, however, that for the wings of equal root beading 
stress, those with greater taper have higher aspect ratios, and hence the 
drag behavior of these wings is, in effect, due to aspect— ratio variations. 

Variation of wing wave— drag coefficient with sweepback for taper 
ratio 0.5 and aspect ratio 2 at a Mach number of 3 (supersonic maximum- 
thickness line) is shown in figure 8. Increased sweepback increases the 
drag rather than decreases the drag as is the case at lower Mach numbers, 
that is, for subsonic maximum-thickness lines. Calculations at other 
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Mack numbers show similar results for the effect of sweepback. It should 
he home in mind, however., that a constant aspect ratio -was maintained; 
that is, increased sweepback was obtained by sliding each section rearward 
rather than rotating the wing panels rearward. 


CONCLUSIONS 


The theoretical investigation of supersonic wave drag of sweptback 
tapered wings at zero lift with thin symmetrical double-^wedge sections 
with maximum thickness at 50 percent chord (rhombic profiles) has been 
extended to include "supersonic" maximum-thickness lines; that Is, the 
flight velocity component normal to the line of maximum thickness is 
supersonic. This condition exists at Mach numbers for which the Mach 
lines have angles of sweep greater than that of the line of maximum thick- 
ness. For purposes of completeness, the results obtained in NACA TN No. Ikk8 
for "subsonic" inaximum-thlckness lines are included in the following 
conclusions . 

1. For wings of constant aspect ratio and given sweepback, taper 
increases the wing wave— drag coefficient at Mach numbers for which, the 
xnaximum-thickness line is substantially subsonic, decreases the drag 
coefficient" in the near— sonic through moderately supersonic range, and 

has negligible effect when the maximun^-thickness line is highly supersonic. 

2. For given sweep and taper ratio, higher aspect ratios reduce the 
wing wave-drag coefficient at Mach numbers for which-the line of maximum 
thickness is substantially subsonic, increase the drag coefficient in 
the intermediate range, and have negligible effect when the maximum— 
thickness line is highly supersonic. 

3. For wings of equal root bending stress and given sweepback, taper 
reduces the wing wave-drag coefficient at Mach numbers for which the 
maximum-thickness line is substantially subsonic, increases the drag 
coefficient in the intermediate range, and has negligible effect when the 
maximunt-thickness line is highly supersonic. 

k. For given taper ratio and aspect ratio, increased sweepback 
reduces the wing waverklrag coefficient at speeds corresponding to subsonic 
maxi mum-thickness lines and increases the drag coefficient when the 
maximumr-thicknees line is supersonic . 

Langley Aeronautical laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., February 25, 19k8 
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APPENDIX A 


EVALUATION OF EQUATIONS (4) FOB SECTION WAVE 
DEAD- AND WING- DRAG EXCLUSIVE OF 

TIP EFFECTS 0 < X, g l) 
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where 

2H1Q 

A = i/nv - 1 
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and 
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m l / 


Wing Drag 
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-1 1 


■— , M — - — — cos ■*- ~ 

- mo) (m-L + j ? 2 m Q 2 - 1 “cr 


(mg - mjjelmQ - 2^0 g 2 , (l - n^nigP^dmQ + 2amg 

+ — ar* COS" — — 

2mg(m 2 - m-J ]/ ^ 2 m 1 2 - 1 0| (“2 " 


(mg + mOdn^ + 2am 1 mg 2 ' (L + m^mgP^dmQ + 2am2 

4=£: — — — _ — =L cos” 

2mg(mg + m^) \JJ 2 ^ - 1 P | ( m l + + 2am l m 2 


vs. 


-1 1 


2 2 
s lEgin^ 


; — -- COS" + — - CO£ 

(ffiQ + Htj) j/e^ 2 - 1 m l p (niQ + mg) [/ p 2 mg 2 - 1 


ka. 2 ia.^pm^ 


(mg 2 - mx 2 )^ P 2 mg 2 - 1 


COS 


.1 2 °“2 


G~ 

ts 

- <±m 


(? V - 0 


Samg 2 ? 


D = 


rtm 0 


\l P%) 2 " 


a 2 m^ d (ding - dniQ + 2amg) 




m l _ “0 


El 


+ 



(2 2\ ,/~2 2 
\™2 - ™0 / (/ P “2 


1 


-1 _1_ 
nigP 
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E = 


*“0 


(/ - 


A . 2 


“0 " m l 


L 1 

+ d^im^ - dmQ + 2 am-^ 


2m 2 m 1 m 0 ^ 

J \l (3 2 m 0 2 - 1 (m-L 2 - Mq 2 ) 


and 


F = 


rr 


jtniQ 


\kW - 


i- 1/ / 

* *■. ✓v/ -i- _ 

a 2 m 0 m 2 dmQ^pani^ - daa g + dm^ 

3^2 *" ^0 ^2 


i2m iV 


Au- 2 


(3™! - “o) - 


-1 1 . a “l m 0 -1 1 

cos ■ — + , — — ■ cos - — 

Pm2 (mQ + \j - 1 ^ m l 


4 aV - 2 


'“l 


COS 


-1 1 


(3 m i - ^ (mo + mi) j/ ? 2 m 0 2 - 1 P “° 
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APPENDIX B 


EVALUATION OF EQUATION ( 6 ) FOE TIP EFFECTS 0 £ A ^ l) 


For aspect-ratio limitations, see equation (5)* 


Section Drag Increment- 


For — < p ^ , 

“O m l 


rtc*,. c 


' dti ^ g - A + B + C 

dm Q 2 (i + j3au) - am^Q 

where A is evaluated in the region — 1 tL- — — % y 'a a™ 

and is equal to 


y(mp - m 1 j - _ ± ym-^l - m^P 2 ) - dm^ 2 (l - + anL^mQ 

m^\jl - m ^ 2 K| 7 5 "o - m l) - “V 1 ! 


- 2 


(7 - am 0 _) 


cosh' 


^ jm^ + am^Q - dm^* 

^vbC^o - y) 


dmQUip^i + Pm-,) - Pam^m 

B is evaluated in the region — — < v < 


m^l + Pmg) 


o y 5 dmQ 


and is equal to 


2 jjmQ - m-^y - amQiiiJ _ x 7 (2m! - “q ' + 2am 1 m p ~ ^q 2 ^- - m i^^) 

iDq \jl - m-L 2 ? 2 2pm 1 |y(in 0 - - am^ 


•, y^Pm.-, - ihq) + ^feaui-, - 
+ (y - dmj cosh "- 1 \ x Z 2 L -Hi ± — 

pm i m o(^o - y ) 
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and C ie evaluated in the region 


dmg (l + - a mg 


1 + 




^ J S dn^ 


and la equal to 

2(5 ( m 2 - ™ o ) - ^.3, ' “oV!) “ *“2^ “ + 

mg /p 2 ^ 2 - 1 P |y^2 - “q) - am 0 m 2j 


-2(7-61^) 


cosh 


-1 


y + m 2 (a - d) 

Pmg (dniQ - j) 


For B > — , use the same formulas as for — < B < JL: hut 
m l “0 ” m l 

change \j 1 - m-^p 2 to \J $ 2 m^ 2 - 1 and cosh - '*' to cos - '*' In the first 
term, only of hoth A and B. 


Wing Drag Increment 


The Increment In -wing wave drag caused h7 the tip is identical^ 
equal to zero for all cases satisf7lng these aspect- ratio limitations. 
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(a) Aspect ratio, 2. 



(b) Aspect ratio, 3. 


Figure 6.- Variation oi wing wave-drag coefficient with Mach number for wings of constant 
aspect ratio. Sweepbsck of maximum-thickness line, 60°. Mach lines are parallel to 
the maximum-thickness line at Mach number 2.0. 
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Straight wmg of infinite aspect ratio 


A = 3.16 




Figure 6.- Variation.’ of wing wave-drag coefficient with Mach number for wings of equal root 
bending stress. Sweep back of marl mum -thickness line, 60°. Mach lines are parallel 
to the marlmnm-thirilniess line at Mach number S.O. 
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Figure 7.- Variation of wing wave-drag coefficient with Mach number for wings of taper ratio 0. 
Sweepback of maximum-thickness line, 60°: Mach lines are parallel to the maximum-thickness 
line at Mach number 2.0. 









Figure 8.- Variation of wing wave -drag coefficient with sweepback of 
maximum -thickness line. Taper ratio, 0.5; aspect ratio, 2; 

Mach number. 3. 





